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Abstract. We consider spatially averaged inhomogeneous universe models and argue 
bJQ, that, already in the absence of sources, an effective scalar field arises through foliating 

and spatially averaging inhomogeneous geometrical curvature invariants of the Einstein 
£Nj ■ vacuum. This scalar field (the 'morphon') acts as an inflaton, if we prescribe a potential 

of some generic form. We show that, for any initially negative average spatial curvature, 
the morphon is driven through an inflationary phase and leads - on average - to a 
ly-j . spatially flat, homogeneous and isotropic universe model, providing initial conditions 

for prc-heating and, by the same mechanism, a possibly natural self-exit. 



1. Introduction 

For about 30 years now, inflation [1] is the main paradigm used to explain the various 
caveats of the Hot Big Bang scenario. In particular, the flatness, smoothness and 
horizon issues [2] are often advocated to be solved by a long enough period of de Sitter- 
like expansion, provoked by of a slow-rolling scalar (multi-)field that dominates over 
the other components of the energy budget of the Universe [3]. The success of the 
concept of inflation stems not only from its simplicity but also from the consequences 
this framework bears, i.e. conservative predictions concerning scale invariance, the 
tensor-to-scalar ratio, Gaussianity of fluctuations, etc. However, as many paradigms 
praised for their consequences, inflation lacks a physical cause. Indeed, despite several 
attempts to justify the global predominance of a yet unobserved fundamental scalar 
field, the existence of the inflaton remains a conundrum to theorists. 

Confronted with this issue, physicists might roughly choose between three distinct 
philosophies. Either a non-standard part of the particle physics spectrum, for instance, 
the Higgs boson [3], or a pseudo Nambu-Goldstone boson [5], could be candidates for 
the fundamental field whose potential part dominated the Early Universe. Or, a known 
part of the physical fields could emulate an effective scalar field via a known process like 
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boson or fermion condensation |JJ E]. Or, a neglected part in the underlying theories, 
due to their special nature or just due to idealizing assumptions, prevented us so far to 
take into account physical effects that could lead to an inflationary era. The idea we 
present here lies within this last line of thought by applying to inflation what has been 
proposed previously for a conservative explanation of the dark energy problem through 
inhomogeneities |S]. In concrete terms we shall provide a physical cause to inflation by 
identifying the inflaton with an effective classical field that describes the averaged bulk 
effect on the dynamics of inhomogeneities (also dubbed backreaction) prior to, during 
and eventually after the inflationary era. 

It is important to stress that the model expounded here pertains to the dynamics 
of backreaction - described by a conservative set of equations and a minimal number 
of physical assumptions - and not to a suitably chosen set of initial conditions. 
Furthermore, as a natural starting point, the inhomogeneities we shall describe hereby 
are due to geometry and a priori not related to the baryonic/dark matter, dark energy 
or radiation content of the Universe [9]; we shall confine ourselves to inhomogeneities 
in the gravitational field only. In other words, this model will demonstrate in simple 
terms how inflation can naturally emerge out of vacuum inhomogeneities. This last term 
does not refer to quantum vacuum fluctuations on a FLRW background (as practically 
always used in the cosmological context), but to the average of classical fluctuations of 
the gravitational field only. 

2. Inhomogeneous universe models 

Einstein's equations are assumed to hold for a four-dimensional tube of the Early 
Universe, featuring the 4— Ricci tensor R^ v = C0- We suggest further to foliate this 
tube into three-dimensional space-like hypersurfaces according to the ADM formalism. 
We here choose a comoving-synchronous line-element (in forthcoming papers the 
embedding issue will be thoroughly addressed), ds 2 = —dt 2 + g^dX 1 dX\ where the 
proper time t labels the hypersurfaces and X % are Gaussian normal coordinates; gij are 
the components of the full inhomogeneous 3— metric of the hypersurfaces of constant 
proper time. In this metric the components of the 4— Ricci tensor can be expressed 
geometrically through the three-dimensional extrinsic curvature of the embedding 
into spacetime, and instrinsic curvature IZij of the hypersurfaces at constant t with the 
following well-known Gaufi-Codazzi-Mainardi relations: 

R 00 = K-K ij K ij , R 0i = K li -K k i]lk , R ij = K ij -K ij -2K ik K k j + KK ij , (1) 

where a vertical slash denotes partial spatial derivative with respect to X 1 , a double 
vertical slash covariant spatial differentiation with respect to the 3— metric, and an 
overdot the covariant time derivative. We end up with the following set of equations 
(supplemented by the defining equation for Ky): 

flu -12 A',,, - A' • K ij /\',, = ; K lt -K% k = ; —k ij +KK i j—2K ik K k j = 7^.(2) 
$ Greek indices run through 0, 1, 2, 3, while latin indices run through 1, 2, 3, and we set c = 8nG = 1. 
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The second equation is the vacuum version of Raychaudhuri's equation, the third set 
of equations are the momentum constraints. Forming the trace of the last equation 
and inserting it into Raychaudhuri's equation, we obtain the vacuum version of the 
Hamiltonian constraint: 1Z + K 2 — K % - K\ = 0, where TZ = TZ\ is the scalar curvature 
with respect to the 3— metric, and K := K\ = — can be interpreted as (minus) the 
local expansion rate of the hypersurfaces. We further propose to spatially average the 
scalar parts of the above equations, defined for any scalar function ty(t, X 1 ) as 

^>{fX l )) v := ±- [ Jdet( 9ij )d 3 X , (3) 

where the volume of an arbitrary compact domain is Vx>(t) := L, ^det(gij)d 3 X. 
Defining a volume scale factor by a v (t) := (Vx>(t)/Vx>(tj)) , and averaging 
Raychaudhuri's equation and the Hamiltonian constraint using the non-commutativity 
relation (true for any scalar S), (S)'x> — (S)x> = (@ S) v — (Q) v (S) v , we obtain the 
following well-known equations [TO] : 

; ^-|-> + & >, (4) 

where H v denotes the (domain-dependent) Hubble rate H v = dx>/ax> = — 1/3 (K) v , 
and A/z?i a (domain-dependent) constant of integration. The kinematical backreaction 
Qd is composed of extrinsic curvature invariants, while is an intrinsic curvature 
invariant and describes the deviation from constant-curvature, 

Q v := (K 2 - K)K\) V - | (K) 2 V ; W v := (TZ) v - ^ . (5) 

For a homogeneous domain the above backreaction terms, being covariantly defined with 
respect to a given spatial embedding, are zero. Therefore they encode the departure 
from spatial homogeneity. The integrability condition connecting the two Eqs.(jl]), reads: 

a- 2 (a 2 v W v ) = -a-\a%Q v ) , (6) 

expressing a combined conservation law for intrinsic curvature and extrinsic fluctuations. 



3. Effective scalar field: the morphon 

We rewrite the set of spatially averaged equations ( H1I6T) and cast it into a Friedmannian 
form with a purely geometrical effective energy-momentum tensor [HJ [12]: 

d v gf s + 3pgf 2 3fc Pi v v v s 

3 — , 6ti v H — - g eS , g cS + iti v ^ eff + p cS ) - U , (7) 

uj) £ (I'd 

where the effective densities are "morphed" by a (minimally coupled) scalar field, the 
morphon $p, thoroughly discussed in [T3] and here defined through: 



QcS-= g = ^v + Uv ; p * := 6 =^v-Ut>, (8) 
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where £ = +/ — 1 for a standard/phantom scalar field, and Uv = Uv{^v) is the self- 
interaction potential. From the above equations we obtain the following correspondence: 

Qv = Ux> — &l ; Wv = -3U V . (9) 

We appreciate that the deviation of the averaged scalar curvature from a constant- 
curvature model Wd is directly proportional to the potential energy; an expanding (or 
contracting) classical vacuum with on average negative scalar curvature deviation (a 
positive potential Uv) can be attributed to a negative potential energy of a morphon 
field. The homogeneous case, Qv = 0, corresponds to a virial equilibrium of the scalar 
field energies. Inserting (Q into the integrability condition fl6]) implies that $x>, for 
$D 7^ 0, obeys the (scale-dependent) Klein-Gordon equation: 

Qv + SHv&v + t-J^-Uv = . (10) 

Averaged universe models obeying the previous set of equations follow, thus, 
Friedmannian kinematics with no fundamental sources, but with an effective scalar 
field perfect fluid that reflects the shape of spatial hypersurfaces and their embedding 
into spacetime. Given a choice of £ and of the potential, the evolution of these models 
is fixed (the governing equations are closed). 

In order to model inflation we impose the constraints (Uv > 0,£ = +1) =>- (Wp < 
0, Wv + 3<2x> < 0): the morphon is canonical and its potential competes with its 
kinetic energy. These conditions still allow for any sign of the kinematical backreaction 
term; rewriting the extrinsic curvature invariants in (jSJ) in terms of the local expansion 
and shear [8], we get that geometries dominated by their expansion fluctuations have 
Qv > -v^> Uv > $!> (hence dv > 0, see fllj), those dominated by their shear 
fluctuations, Qv < 0; and homogeneous spacetimes obey Qv = on all scales. 



4. A morphonic inflaton 

We offer the idea that the morphon can formally play the role of the inflaton, and that 
(unavoidable) curvature inhomogeneities occurring at one point of the Universe's history 
could be the actual cause of a de Sitter-like era prior to the baryon nucleosynthesis 
era, provided they are governed by an appropriate potential. In the average quantities 
dictionary, inflation can be read off in terms of the effective first slow-roll parameter 

a v > e:= -Hv/H v < 1 . (11) 

Some remarks are already necessary at this stage in order to assess our paradigm and to 
distinguish it from the usual "FLRW+fundamental sources" inflationary picture. First, 
since Eqs. (Hj) only functionally depend on a metric that needs not to be specified, 
requiring dv > does not imply that the comoving Hubble distance 1/aH decreases; 
however, such is the essence of the fitting problem to find the best FLRW fitting model 
to a lumpy universe [HI EES]- Secondly, a fair solution of the flatness and smoothness 
problems would necessitate to implement averaging on the light cone, or would require 
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an explicit inhomogeneous metric. Lacking the latter in the present context, we claim 
that realizing e < 1 during a sufficient number of e-folds is the best condition to address 
the dynamics of inflation. Third, so as to legitimate the neglect of any averaged matter, 
we refer to [26], where it is shown with the same class of assumptions (a flat enough 
curvature invariant W-p, see (112j) below) that the corresponding cosmological parameter 

= (p)v/3Ht> dives under a percent within a few H^, l {ti) times. Therefore a large 
enough initial "Hubble parameter" Hx>{ti) consigns the matter content to oblivion. 

A suitable potential Ut> has to be chosen in order to identify the morphon as a 
trustworthy inflaton. This identification implies that the cause (the inhomogeneities) 
and the consequence (the smoothness) of the present model are at odds. Therefore, 
a meaningful potential needs not only to describe how inhomogeneities can induce 
inflation, but also how they can disappear throughout this process. Translated to 
the standard inflationary paradigm, this condition can be easily implemented by any 
potential that possesses a minimum downwards which the inflaton rolls. One of the 
simplest examples, which has been extensively studied in the context of chaotic inflation 
[TT] . is a potential of the Ginzburg-Landau form: 

U% L = U ($,-<f$f/*t. (12) 

This quartic potential can also be related to a fundamental Higgs field. However, even 
if such a scalar field is fundamental, there is the possibility that no Higgs particle is 
involved - as in our case - e.g. [18]. Contrary to standard inflation one must here raise 
the issue of the "reality" of such a potential, since we postulate that U-p is actually 
given by the averaged scalar curvature ([9]). More precisely, one should question the 
assumption of a slow-roll period, which corresponds to a nearly constant value of the 
averaged curvature. At this stage we cannot provide a proof for this possibility, e.g. 
through a suitable exact solution or through a well-defined approximation for the local 
evolution of curvature, but we can summarize plausible hints that the actual physical 
properties of an intrinsic curvature distribution could comply with the expected: first, 
there are indications from perturbation theory that the perturbative expansion is led by 
a term Wp oc a^ 1 , when putting the Friedmannian zero-order curvature to zero (so that 
the scenario is exclusively governed by backreaction, i.e. the curvature-deviation from 
a flat model), and this perturbative expansion can be extended up to including even 
a constant term [19]; secondly, a concrete modeling of the non-perturbative aspects of 
the curvature distribution in a multiscale analysis reveals that, even if subdomains are 
evolving according to the term oc a^ 1 , the variance between subdomains can lead 
to a de Sitter-like behavior on large scales [16] . These results, though not being a proof 
of the actual choice we made (say, we employ (fl2"|) as an illustrative example), open the 
possibility that there exist some configurations of the geometrical inhomogeneities that 
enable a slowly varying curvature, that is a sufficiently flat potential. Obviously, we here 
are in the same situation as for an explanation of the dark energy phenomenon through 
inhomogeneities, an ongoing research field that is qualitatively well-understood as for 
the mechanism but is not yet quantitatively conclusive. 
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Figure 1. The Ginzburg-Landau potential in arbitrary units and the possible initial 
conditions as well as their physical meaning. All conditions possess some curvature 
Wv t < 0. The arrows schematically indicate the amplitude of the morphon's initial 
speed $£>.. In the order of the points (from left to right): the first two points dominated 
by shear fluctuations (red, green) are obtained for Q-p. < <£>|,. > 2(iJ|,. +&£>.); the 
next points dominated by expansion fluctuations (pink, blue) for < 2(_ff|,. + fcx>. ), 
where the de Sitter-A equivalent case has a stiff morphon $£>. = 0; the homogeneous 
case (last point, orange) is obtained for = 2(U|,. + fcx>,). 



Once the minimum $ is fixed, the evolution of the morphon, given the integrability 
condition fl6]), is practically independent of the initial conditions^ [20]. In Fig. [1] we 
show how all acceptable initial conditions are reinterpreted in terms of the curvature, 
and expansion/shear fluctuations. For any of these types we find the behavior shown 
in the inset of Fig. [2j instead of simple inflation e < 1, the choice of U% L gives us 
slow-roll inflation t C 1 during which the energy density budget is dominated by 
the curvature term Vt^ = — >Vr>/6if|, — 3/2, and a negative kinematical backreaction 
density fi^ = -Qv/QH^ ~ -1/2, i.e. by expansion variance, whatever the initial 
value of the homogeneous part of the curvature, Q^ 1 = —kj ) Ja\ ) R\ ) = 1 — Qq — Q^, 
is. At the end of the inflationary era, the morphon oscillates down its potential and, 
because of the order of the polynomial (fl2|) . the effective scale factor behaves as if it 
were dominated by dust matter, a-p oc t 2 ^ 3 , and the effective deceleration parameter 
q° = —a v ,{ax>H'^ ) )~ 1 = 2VLq oscillates around q v = 1/2. While oscillates around 
3/4 and VLq around 1/4, and Qd tend to vanish, i.e. our model generically tends 

§ Uq determines the inflation's onset time and influences also its duration. Due to the necessity 
to recover the CMB temperature fluctuations, "classical" inflationary models suffer from fine-tuning 
issues, once perturbation theory is performed on the inflaton. Such is not the case here, since our model 
tackles the bulk effect of inflation within a highly non-linear regime, without invoking perturbation 
theory Hence, though a lack, the absence of perturbative constraints allows us to get rid of fine-tuning. 
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Figure 2. The energy densities fig (solid, red), Ciy^ (dashed, blue), and flf 
(dotdashed, green) for one of the cases of Fig. Q] as a function of the number of e-folds 
N = In a-p . The initial value of the homogeneous part of the intrinsic curvature has 
been taken arbitrary large (e.g. QjP' ~ 1) to underline the fact that it vanishes anyway 
after a few e-folds. In the inset, the slow-roll parameter e for the same configuration 
is shown. Since e = f2j[? + 4>|>/2 always holds, e ~ 1 in the pre-inflation era when 
ttf > and e ~ during inflation when U% L > $f,/2. 



to be - on average - homogeneous and quasi-isotropic. 

Let us sum up the evolution scenario driven by kinematical backreaction Qx>: 
thanks to a favorable configuration of the latter and of the average curvature, a) a 
period of accelerated expansion kicks out (Q-p rises more rapidly than H%,), b) an 
inflationary era follows during which the expansion variance stays positive and almost 
constant (Qx> — 3iJ|, ~ const.), c) the outcome is a globally inflated universe, hence 
exponentially smooth, where gravity and shear pull back, causing on time-average a 
negative kinematical backreaction Qx> — —2/3t 2 . Hence, the inhomogeneities wash out 
at the end of the process, so does the acceleration, hereby providing a natural graceful 
exit, univocally based, should we stress, on mere gravitation. 

5. Discussion 

Looking at a portion of the classical vacuum within Einstein's theory, we obtained, 
by foliating spacetime and spatially averaging the scalar parts of the Einstein vacuum 
equations, a Friedmann-like evolution of the scale factor, however, driven by an effective 
scalar field. Starting with a sufficiently flat potential, the backreaction mechanism 
initiates a metamorphosis of the geometrical properties of space that goes along with 
an exponential inflationary phase. This dynamics tends to flatten the averaged scalar 
curvature and suppresses fluctuations. Although the present scenario just prescribes the 
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initial conditions for a follow-up pre-heating, the same mechanism sets the conditions 
to provoke an intrinsic exit scenario. Note also that, while acceleration was associated 
to the formation of inhomogeneities in the Late Universe, e.g. [21], the application of 
the same mechanism to the Early Universe tends to homogenize spatial hypersurfaces. 

Let us list some of the immediate consequences of this mechanism. First, due to 
its classical nature, the inflaton's mass is no longer limited by the Planck/SUSY scale 
(see e.g. [22] ). Secondly, inflation is already possible for the unmodified Hilbert action, 
which of course does not exclude the need for improvements of Einstein's theory. Third, 
inflation, often dubbed to be unsustainable in inhomogeneous spacetimes [23] (though 
see e.g. [21]), could occur despite the natural chaotic (inhomogeneous) initial conditions 
invoked by the theory [25] (the reason for this behavior is essentially that we do not 
study the stability of the fluctuations on a homogeneous reference background but that 
of a background-free general average). 

Once this average correctly computed, one can address the issue of fluctuations 
which one expects to be the seeds of large-scale structure and of the CMB spectrum; 
the absence of such a prediction is certainly a dearth that we should overcome (by 
formulating of a fluctuation theory about the average). In a joined paper [26], we 
improved the present model by proving that a foliation into flat space sections is unstable 
and is attracted by an inhomogeneous negative-curvature state, and also by proving that 
the "empty" model we considered is a generic attractor of "filled" models under some 
conditions. In a forthcoming paper, we shall address a more general model that includes 
matter, radiation and fundamental scalar field inhomogeneities. If sources are present, 
there are interesting interactions with a morphon field, the latter being always present 
in the case of an inhomogeneous cosmology. 

The presented scenario points to a huge potential of studying scalar field models 
in the Early Universe. We gave the simplest conceivable model that generates inflation 
out of the classical inhomogeneous vacuum. Contrary to the quasi-Newtonian standard 
picture, we render the bulk effect of curvature responsible for inflation. In other words, 
more than being acquitted of preventing inflation, inhomogeneities, when treated non- 
perturbatively, could be the actual cause of it. 
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